Biophysical Journal

BACKGROUND In line with International Federation of Clinical Chemistry and Laboratory Medicine (IFCC) recommendations for high-sensitivity cardiac troponin assays, analytical validation and reference limit assessments are required to confirm that an assay meets performance criteria. This study evaluated the analytical performance and established the 99th percentile upper reference limit (URL) for the SPINCHIP High-Sensitivity Cardiac Troponin I (SPINCHIP hs-cTnI) point-of-care assay. METHODS Analytical performance characteristics, including the limit of blank (LoB), limit of detection (LoD), and limit of quantification (LoQ), were assessed. Additionally, 1,053 plasma samples and 1,055 whole-blood samples were used to determine the URL. Imprecision around the 99th percentile URL was evaluated as part of the analytical validation. High-sensitivity criteria were assessed by confirming measurable cTnI in [&ge;]50% of healthy individuals (n=432 plasma; n=431 whole blood) and achieving imprecision <10% at the 99th percentile (plasma, n=960; whole blood, n=480). RESULTS SPINCHIP hs-cTnI demonstrated a LoB of 0.3 ng/L; LoDs of 0.8 ng/L (plasma) and 0.9 ng/L (whole blood); and LoQs of 1.1 ng/L (plasma) and 1.4 ng/L (whole blood). The analytical measuring range was 1.1-9,000 ng/L. Imprecision at the common 99th percentile URL (14 ng/L) was 5.8%; for men (URL=16 ng/L) 5.6% and for women (URL=10 ng/L) 6.3%. Greater than 85.2% (94.0% and 76.1% in men and women, respectively) of healthy individuals showed measurable cTnI above the LoD. CONCLUSIONS The SPINCHIP hs-cTnI assay meets the IFCC high-sensitivity requirements, demonstrating <10% imprecision at the 99th percentile, reliable low-concentration precision and cTnI detection in more than half of healthy individuals.

Background Hypertension is a major modifiable risk factor for atrial fibrillation (AF), yet blood pressure (BP) control remains suboptimal in older U.S. adults. Objectives This study evaluated how improve systolic BP (SBP) control could affect incident AF, downstream AF ablation demand, Medicare savings, and hospital revenue. Methods A population-based modelling framework was developed to estimate mortality and incident AF hazards across SBP strata: <120, 120-139, 140-159, and ?160 mm/Hg. AF incidence in the SBP <120 mmHg group was set at 2.2 per 1,000 person-year, with hazard ratios of 1.17, 1.42 and 1.64 applied to higher SBP strata. We assumed 25% of incident AF patients would undergo ablation, with a 7.2% complication rate. AF prevalence was projected to increase by 4.6% annually over 10 years. Medicare savings and hospital revenue foregone were estimated under varying procedure cost and contribution-margin assumptions. Results Higher SBP was associated with greater hazards of death and incident AF. Improved SBP control reduced projected AF incidence and ablation demand. Over 10 years, cumulative Medicare savings were projected at $8.7B-$10.9B across the full modelled population. However, reduced ablation volume translated into hospital revenue foregone, ranging from $75M to $377M in the first year, and approximately $1.03B-$5.2B cumulatively over 10 years. Conclusions Improved SBP control may reduce AF incidence, prevent avoidable invasive ablation procedures, relieve pressure on surgical waitlists, and generate substantial Medicare savings. However, these benefits may reduce hospital procedural revenue, highlighting a misalignment between prevention-oriented care and fee-for-service reimbursement incentives.

Minimally invasive endovascular procedures offer reduced surgical trauma, shorter recovery times, and improved outcomes, but rely on 2D fluoroscopic X-ray imaging, which provides limited depth perception and exposes patients and clinicians to ionizing radiation. Here we present an augmented reality (AR) system that fuses intravascular ultrasound (IVUS) and electromagnetic (EM) position tracking with preoperative computed tomography (CT) to produce an anatomically accurate, deformation-corrected navigational reference. A robotic device performs ECG-gated pullback of the IVUS probe, capturing 4D aortic motion across the cardiac cycle. We introduce a deep learning architecture for extracting vascular lumen boundaries and side-branch orifices from artifact-prone IVUS streams, and a semantically driven non-rigid CT-IVUS fusion pipeline robust to false positive landmarks. We evaluate the platform with trained surgeons in benchtop phantom studies and in-vivo ovine models, and demonstrate its application to fenestrated endovascular aneurysm repair (FEVAR). Compared to fluoroscopy alone, AR guidance significantly reduces cannulation time, radiation exposure, and cognitive workload, while improving procedural efficiency and safety. Our IVUS-EM and CT aortic datasets are released open source.

ABSTRACT Background: Athletes who return to soccer after anterior cruciate ligament reconstruction (ACLR) remain at elevated risk of secondary injury despite meeting conventional discharge criteria, and neuromuscular deficits in the reconstructed limb are known to be exposed by fatigue. Objective: To determine whether match-play fatigue differentially affects muscle stiffness, countermovement jump (CMJ) force symmetry, and rate of force development (RFD) asymmetry between soccer players with a history of ACLR and uninjured teammates. Methods: A prospective, cross-sectional, matched-control study enrolled 128 competitive soccer players (64 ACLR, 6-22 months post-surgery; 64 uninjured controls) across five recruitment waves (February-June 2026). Bilateral CMJ peak vertical force, jump height, RFD, and myotonometric stiffness of the rectus femoris (RF), vastus medialis (VM), and biceps femoris (BF) were recorded immediately before and after a standardized competitive match. Fatigue was quantified from second-half heart rate (percentage of age-predicted maximum) and end-match rating of perceived exertion (RPE). Within-group pre-to-post changes were evaluated with paired t-tests, between-group differences in the magnitude of change with independent-samples t-tests, and associations between fatigue indices and asymmetry changes with Pearson correlations. Results: Match play reduced CMJ limb symmetry index (LSI) in both groups, but the decline was more than three-fold greater in the ACLR group, 92.6% (SD 5.4%) to 85.1% (SD 7.1%), than in control group, 97.3% (SD 3.9%) to 95.0% (SD 4.2%), group-by-time difference, p < 0.001, (d = 0.64). RFD asymmetry approximately doubled in the ACLR group, 10.6% (SD 4.1%) to 17.6% (SD 6.5%), compared with a smaller rise in control group, 4.6% (SD 2.4%) to 6.3% (SD 3.7%); p < 0.001, d = 0.77). Involved-limb stiffness losses in the ACLR group exceeded those of controls for the RF (-21.2 vs. -9.2 N/m, p < 0.001), VM (-17.7 vs. -6.1 N/m, p < 0.001), and BF (-13.3 vs. -6.6 N/m, p < 0.001), whereas uninvolved-limb stiffness losses did not differ between groups (all p > 0.05). Fatigue markers (heart rate, RPE) were not significantly correlated with the magnitude of individual asymmetry change (|r| [&le;] 0.18, p > 0.15). Conclusions: In competitive soccer players 6-22 months after ACLR, match-play fatigue selectively compromises stiffness and explosive force output of the reconstructed limb, widening inter-limb asymmetries beyond what is seen in uninjured teammates, even though global cardiovascular and perceptual fatigue were comparable between groups. These findings suggest that return-to-sport testing performed only in a rested state may underestimate residual neuromuscular deficits, and support fatigue-inclusive assessment protocols before athletes are cleared for unrestricted competition. Abbreviations: ACL: anterior cruciate ligament, ACLR: anterior cruciate ligament reconstruction, BF: biceps femoris, CMJ: countermovement jump, HRmax: maximum heart rate, LSI: limb symmetry index, RF: rectus femoris, RFD: rate of force development, RPE: rating of perceived exertion, RTS: return to sport, VM: vastus medialis, SD: standard deviation. Keywords: Anterior cruciate ligament reconstruction, muscle fatigue, muscle stiffness, countermovement jump, limb symmetry index, rate of force development, soccer, return to sport.

Background: North Carolina Medicaid implemented the Pregnancy Medical Home program to improve access to high-quality maternity care and reduce the risk of adverse perinatal outcomes. Program recipients receive a prenatal risk screening form, originally intended to identify those at high risk of preterm birth and low birth weight, that includes an assessment of social and clinical factors. While prior studies have evaluated whether risk screening can identify pregnancies with higher risk of adverse neonatal outcomes, less is known about the relationship between programmatic risk-stratification and adverse maternal outcomes. Objective: To assess the use of a prenatal risk screen among pregnant Medicaid beneficiaries to identify those at risk of an adverse maternal event. Study design: Linked Medicaid hospital claims, live birth records, and risk screen data from the Pregnancy Medical Home program were used to identify risk factors for adverse maternal events among individuals who gave birth to a liveborn infant in North Carolina between 2014 and 2019. Only those with completed risk screens (75%) were included in the analysis. We used random forest classification to select variables for a multivariable prediction model. We used Poisson regression to model the association between adverse maternal events and selected demographic, psychosocial, clinical, and historical pregnancy characteristics. Adverse maternal events occurring at birth and up to six weeks postpartum included severe maternal morbidity, maternal intensive care unit admission, prolonged birth hospitalization, and postpartum readmissions. Results: A total of 205,916 births met inclusion criteria for this analysis. During the study period, 3.0% of Medicaid beneficiaries had an adverse maternal event occurring between birth and up to six weeks postpartum, including, 0.6% with severe maternal morbidity, 0.9% with an intensive care unit admission at birth, and 1.5% with a prolonged birth hospitalization or postpartum readmission. Maternal age greater than 25 years, Black race, being overweight or obese, smoking, chronic diseases (diabetes, hypertension, mental illness), and pregnancy history characteristics (nulliparity, history of preterm birth, history of hypertensive disorders of pregnancy or gestational diabetes) were associated with an increased risk of adverse maternal events. Modeled together, however, risk factors from the risk form were poorly predictive of the composite outcome. The final model had an Area Under the Curve (AUC) of 0.63 with an optimal sensitivity of 56% and specificity of 63%. Conclusion: Care management during pregnancy is an increasingly relevant topic in public health and prenatal care in the United States. The North Carolina Pregnancy Medical Home is a long-standing and robust Medicaid program that can serve as a model for design and implementation. While this program has effectively designed risk-stratification to identify pregnant people at risk of poor neonatal outcomes who benefit from care management, the risk screen poorly identifies pregnant people at risk of adverse maternal outcomes. Care coordination programs are often designed to optimize neonatal outcomes, and this study highlights the need to center and balance maternal health along with neonatal outcomes to address the needs of a very vulnerable population.

Introduction: As healthcare systems increasingly shift toward prevention and community-based care, the demand for physicians in extramural specialties continues to grow. Yet, a mismatch persists between workforce needs and medical students career aspirations. Little is known about how medical students and trainees develop an interest in extramural specialties such as youth health care (YHC). This study explores trainees trajectories toward becoming a youth health care physician (YHCP). Methods: We conducted a qualitative study using semi-structured online interviews with fourteen YHCPs in training. We combined an inductive and deductive approach, applying the person-environment (PE) fit framework to explore participants evolving experiences of fit and misfit. Results: Participants described growing misfit with clinical culture of medical school (i.e. the hidden curriculum), particularly during hospital-based clerkships, combined with limited exposure to YHC. For some, this misfit extended to doubts about becoming a doctor. Over time, participants developed a sense of fit and belonging within YHC, either directly or after exploring other specialties including extramural specialties. Discussion: These findings reframe specialty choice as a longitudinal search for belonging and alignment, in which trainees iteratively explore, evaluate, and refine their sense of fit across contexts. Clerkships serve as key sites for testing fit, yet also expose learners to the clinical culture, including the hidden curriculum. Broadening exposure and supporting reflective fit processes may encourage more medical students to choose extramural specialties, ultimately fostering a more balanced and sustainable alignment of the medical workforce.

As machine learning and artificial intelligence systems are increasingly used in healthcare, rigorous evaluation of their classification performance has become critical. The F1 and F{beta} scores are widely adopted metrics for assessing performance in imbalanced biomedical data. Recently, we introduced psF1, a unified statistical framework for inference and study design for single and comparative F1 and F{beta} scores under the assumption of independent classifiers. In practice, however, benchmarking two classifiers on the same dataset creates a correlated paired setting. Ignoring this intrinsic dependency leads to overestimation of the standard error and a substantial loss of statistical power. To address this, we develop psF1pair, an advanced framework for statistical inference and power analysis that explicitly accounts for correlations between classifier pairs. Extensive simulation studies demonstrate the performance of psF1pair, and its utility is further illustrated through application to a real-world imaging classification system. As expected, higher correlation between classifiers yields narrower confidence intervals and enhanced statistical power. A freely available R package is provided to facilitate implementation, supporting accurate evaluation and study design for predictive and classification models in biomedical research.

Background: Apical sparing of left ventricular longitudinal strain (LS) is an echocardiographic clue to cardiac amyloidosis but may also occur in hypertensive heart disease (HHD). Objectives: To determine whether apical sparing in HHD is associated with regional left ventricular wall stress estimated according to Laplace's law. Methods: We retrospectively studied 1,559 patients with HHD, 47 with light-chain cardiac amyloidosis (ALCA), and 409 normotensive controls. Artificial intelligence-assisted echocardiography quantified segmental LS, wall thickness, and cavity radius at the basal, midventricular, and apical levels. Wall stress was estimated as mean blood pressure (MBP) x radius/(2 x wall thickness). Apical sparing was defined as a relative regional strain ratio (RRSR)[&ge;]1.0. Results: Apical sparing was present in 14 patients with HHD (0.9%), 13 with ALCA (27.7%), and no controls. Among HHD patients with apical sparing, RRSR decreased from 1.11{+/-}0.13 to 0.72{+/-}0.10 after antihypertensive treatment (P<0.001), accompanied by reduced wall stress and improved basal and midventricular LS, with resolution of apical sparing in all 14 patients. In the overall HHD cohort, changes in MBP and left ventricular mass index were independently associated with changes in RRSR. In an exploratory analysis of HHD patients with apical sparing, a reduction in basal wall stress was associated with a reduction in RRSR ({beta}=0.267 for {bigtriangleup}RRSRx100, 95% CI 0.023-0.511; P=0.036). In ALCA, favorable hematologic response was the only determinant of RRSR reduction. Conclusions: Apical sparing in HHD was uncommon but reversible and may represent a load-sensitive deformation pattern associated with regional wall stress, consistent with Laplace's law.

CRISPR-Cas12a nucleases are classically activated through CRISPR RNA (crRNA) guided and PAM-dependent target recognition, which together establish a canonical heteroduplex associated with nuclease activation. Here we identify a crRNA- and PAM-independent activation pathway for Cas12a that reveals previously unrecognized conformational plasticity within its nucleic acid recognition interface. We show that short RNAs can directly occupy the canonical crRNA-binding channel and trigger a catalytically competent trans cleavage state in the absence of PAM recognition or canonical R-loop formation. Biochemical assays indicate that short RNAs bind the crRNA-binding channel and are competitively displaced by cognate crRNA, consistent with binding at a conserved nucleic acid-binding interface. Cryo-electron microscopy (cryo-EM) further reveals that Cas12a maintains its global catalytic architecture while exhibiting loss of canonical PAM-dependent stabilization and increased flexibility of the RuvC lid, alongside accommodation of a noncanonical RNA-DNA hybrid with inverted polarity relative to the crRNA-target duplex. This crRNA-independent activation pathway enables programmable, amplification-free detection of DNA and RNA targets independent of canonical guide-mediated recognition. Together, these findings define an alternative activation geometry for Cas12a and expand models of Class 2 CRISPR-Cas effector activation beyond crRNA- and PAM-directed recognition.

The X-linked G-protein coupled receptor GPR174 is highly expressed in T and B lymphocytes and has immunoregulatory roles in mice, but its function in humans is unknown. We describe a cohort of six individuals who have function-disrupting variants in GPR174 and a clinical phenotype of lymphadenopathy and autoimmunity. Histological analysis of two patient lymph nodes revealed necrotizing lymphadenitis and lymphoproliferation resembling Kikuchi-Fujimoto disease. In-depth analysis of three patients and related carriers revealed overaccumulation of CD8 terminally differentiated effector memory cells re-expressing CD45RA (TEMRA). Patient cells and GPR174-deficient CD8 T cells generated from controls showed less repression of proliferation by the GPR174 ligand lysophosphatidylserine (lysoPS) and an effector-biased gene expression program. GPR174-deficient CD4 T cells were resistant to lysoPS-mediated suppression of IL2 production. In mice, chronic viral infection led to over-accumulation of GPR174-deficient effector CD8 T cells. We describe an inborn error of immunity associated with dysregulated lymphocyte responses that we propose predisposes to exaggerated lymphoproliferation and autoimmunity following viral infection.

Background: Monoclonal antibodies have emerged as a prophylactic strategy to prevent symptomatic SARS-CoV-2 infection in immunocompromised individuals. However, the evolutionary and clinical implications of breakthrough infections under this regime remain unclear. Methods: A male in their 80s with a haematological/oncological diagnosis received a 2000 mg intravenous infusion of sotrovimab in March 2023 and was diagnosed with COVID-19 by RT-qPCR from a nasopharyngeal swab in August 2023. Weekly samples (n=24) were collected through February 2024 (171 days). All samples underwent whole-genome sequencing, with select mutations subjected to functional assessment. Findings: Sequencing identified the GE.1 lineage at all timepoints. An intra-host recombination event in ORF1ab (positions 8942-12458) was detected prior to 23 weeks post-detection, followed by a 14-fold increase in viral load (7.42e+06 to 1.00e+08 RNA copies/mL) and a marked shift in the viral population. E340D, a sotrovimab resistance mutation, was detected at low abundance (46%) within the first week post-infection, fluctuated over time, and was nearly fixed by week 15 (107 days) post-detection. We assessed five spike mutations - V36M, S98F, and V213G in the N-terminal domain, Y505P in the receptor-binding domain, and P681Q near the S1/S2 cleavage site - and additionally evaluated the impact of E340D. V36M conferred the highest infectivity across all cell lines, with the most significant effect in low-TMPRSS2 cells. While all mutations showed enhanced infectivity with the addition of E340D, the effect was most pronounced in mutations with lower baseline infectivity. The addition of E340D significantly decreased relative neutralizing titres for V36M, S98F, and V213G, enabling escape from neutralizing antibodies in XBB-responsive individuals, illustrating an enhanced phenotypic advantage. Patient neutralizing activity was absent pre-sotrovimab, and sotrovimab-induced neutralization was further compromised by selection of E340D. Interpretation: Sotrovimab pre-exposure prophylaxis in an immunocompromised patient did not prevent SARS-CoV-2 infection, and selected for resistant mutation E340D, with unexpected fitness consequences across non-receptor binding domain spike regions.

Purpose: Nonresponse to routinely collected patient-reported outcome measures (PROMs) threatens the representativeness of aggregated data. We characterized patient-, provider-, and clinic-level factors associated with PROMIS Global-10 nonresponse in routine radiation oncology care. Methods: In this retrospective cohort study, all adults seen at five Mass General Brigham radiation oncology clinics over one year were included. The primary outcome was patient-level nonresponse, defined as never completing the portal-administered Global-10 versus completing it at least once. Using iterative mixed-effects logistic regression, we modeled patient-, provider-, and clinic-level factors. Results: Among 12,214 patients, 71 providers, and five clinics, patient- and appointment-level response rates were 35.4% and 10.9%, with patient-level response ranging nearly fivefold across clinics (12.8% to 66.2%). In Model 1, male sex, lower education, not working, and recent surgery had higher odds of nonresponse, and longer time since diagnosis lower odds. After provider- and clinic-level factors were added, patient sex, education, and employment became nonsignificant, whereas recent surgery (adjusted odds ratio [aOR] 1.97) and longer time since diagnosis (aOR 0.46 for >12 months) persisted. A provider's historical collection rate was protective but attenuated at the clinic level. There, a later program launch (aOR 0.29) and higher historical collection rate (aOR 0.79) correlated with lower nonresponse, whereas academic versus community setting did not. Conclusions: Nonresponse to routinely collected PROMs is a multilevel phenomenon driven substantially by clinic-level implementation factors, not patient characteristics alone. Because response rate is only a proxy for representativeness, PROMs programs and PRO-based performance measures should prioritize representative collection over volume.

Controlled human infection models (CHIMs) are increasingly used in infectious disease research to study pathogen dynamics and evaluate interventions under controlled conditions. However, these studies are resource-intensive and involve ethical and safety constraints, making efficient study design critical. Dose-finding is a key early component in CHIMs, where the aim is to identify a challenge dose that achieves a target infection probability. Traditional rule-based designs are commonly used but can be inefficient, motivating the use of model-based adaptive approaches such as the Bayesian Continual Reassessment Method (CRM). Although CRM has been extensively studied and widely adopted in Phase I oncology trials for identifying the maximum tolerated dose of therapeutics, its application in CHIM settings remains limited, particularly when the endpoint of interest is infection. This tutorial provides step-by-step guidance for implementing a Bayesian CRM in dose-finding CHIMs, using an oropharyngeal Neisseria gonorrhoeae challenge as a motivating case study. The framework outlines key design components, including dose-grid specification, dose-response model, prior elicitation, Bayesian updating, decision rules, and stopping criteria, with particular emphasis on a clinically interpretable parameterisation. Trial operating characteristics are evaluated through simulation studies under multiple dose-response scenarios and prior-predictive analyses, and compared with a commonly used '3+3' type rule-based design. This work highlights the advantages of Bayesian model-based designs for dose-finding in CHIMs over classic rule-based designs and provides a structured, reproducible framework for implementing CRM, supporting their application in future CHIM studies.

Background and Aims The optimal treatment for first episodes and first recurrences of Clostridioides difficile infections (CDI) is unknown and there is emerging evidence for pulse and taper (P-T) regimens. Therefore, we sought to estimate the relative efficacy of treatment options. Methods MEDLINE and CENTRAL were searched from database inception to May 21, 2025 and unpublished conference abstracts were searched from recent infectious disease conferences. RCTs on the treatment of first episodes or first recurrences of CDI comparing fixed-dose or P-T regimens of fidaxomicin or vancomycin were included. The primary and secondary outcomes were 40- and 56-day CDI recurrence, respectively. A random-effects network meta-analysis on the risk ratio (RR) scale was conducted using a standard regimen (10-14 days) of vancomycin as the comparator. Treatments were ranked using the surface under the cumulative ranking curve (SUCRA). Results 8 RCTs were included comprising a total of 2181 patients. For 40-day recurrence, fidaxomicin P-T had the highest probability of ranking best (RR=0.10, 95%Confidence Interval [95%CI]=0.10-0.49, SUCRA=1.00), followed by vancomycin P-T (RR=0.49, 95%CI=0.32-0.76, SUCRA=0.61), fixed-dose fidaxomicin (RR=0.61, 95%CI=0.49-0.76, SUCRA=0.39), and, finally, fixed-dose of vancomycin (SUCRA=0.00). The treatments ranked in the same order for 56-day recurrence, though only 3 RCTs reported on this timepoint. Conclusion Vancomycin P-T, fidaxomicin P-T, and fixed-dose fidaxomicin were all superior to a fixed-dose vancomycin. Head-to-head comparative effectiveness RCTs are needed to quantify their relative effect sizes of and impact on long-term prevention of recurrent CDI.

Background The 2022-2023 global mpox outbreak highlighted the critical need for robust genomic surveillance capabilities to track mpox virus (MPXV) evolution and transmission dynamics. Methods Building upon our established SARS-CoV-2 sequencing infrastructure, we implemented a Molecular Loop probe-based long-read sequencing approach using Pacific Biosciences Sequel II technology for comprehensive MPXV genomic surveillance across the United States (US). From August 2024 to June 2025, we generated 326 high-quality whole genome sequences from residual mpox-positive clinical specimens collected by Labcorp across all 10 US Department of Health and Human Services regions. Results Our analysis identified two samples containing clade Ib MPXV in January and June 2025 and captured shifting trends in clade IIb diversity, with 13 distinct lineages observed. We also identified multiple instances of large (~1.6-17.6kb) deletions proximal to the inverted terminal repeats in clade IIb genomes. APOBEC3 mutation analysis indicated substantial evidence of human-to-human transmission among both clades. Further, we observed significantly higher APOBEC3-associated SNPs per kilobase (P<0.001) in clade IIb genomic variable regions relative to their central conserved region. Our assay exhibited strong reproducibility across biological replicates from individual patients and accuracy was confirmed via parallel sequencing of select specimens by US Centers for Disease Control and Prevention (CDC) using metagenomic sequencing. We also demonstrated via custom simulation that our assay discriminates all known MPXV clades and lineages, including those we have not observed in the US. Conclusions Our integrated nationwide surveillance system facilitates real-time genomic tracking of outbreak evolution, with demonstrated capacity across SARS-CoV-2 and MPXV, positioning this platform for rapid deployment during future pathogen emergence.

Stochastic epidemic models are a cornerstone of infectious disease epidemiology and are often used to study intervention scenarios. However, large run-to-run variability can make intervention effects difficult to estimate precisely. We revisit the epidemic Sellke construction, which assigns each individual an infection threshold for the cumulative infection hazard such that, conditional on the thresholds, the epidemic trajectory becomes deterministic. This enables coupling of simulations with and without an intervention, yielding low-variance effect estimates even when outcomes such as final size or peak incidence vary widely between runs. We develop an exact, event-driven implementation that maintains infection and recovery events in priority queues. Cumulative infection-hazard updates require O(log N) time per event, yielding overall complexity O(Elog N) for E events in a population of size N. The implementation achieves computational performance comparable to the classical Gillespie algorithm while naturally accommodating non-Markovian infectious periods and complex infectiousness profiles. We illustrate the approach using distance-dependent spread of avian influenza between poultry farms in the Netherlands and a multilayer population with households, schools, and workplaces. In both examples, coupling enables efficient within-run comparisons of intervention scenarios across stochastic realisations.

Background: Opioid use disorder (OUD) is heritable, yet most genome-wide association studies (GWAS) have focused on European populations, leaving the genetic architecture of OUD in non-European populations underexplored. Methods: We conducted GWAS of OUD across three ancestries using electronic health records and genomic data from 52,357 All of Us Research Program participants (8,912 cases; 43,445 matched opioid-exposed controls; 48.5% female). Participants were stratified into European (EUR), African (AFR), and Admixed American (AMR) ancestry groups for logistic regression GWAS, with independent replication in the Million Veteran Program. We then applied the deep-learning model AlphaGenome to predict the tissue-specific transcriptomic and splicing consequences of top risk variants across 13 reward-pathway brain regions. Results: We identified and replicated a novel DDX6 risk locus, alongside established OPRM1 and FURIN signals. AlphaGenome predicted the DDX6 regulatory allele downregulates the stress-resistance gene FOXR1 in the nucleus accumbens, while the protective OPRM1 variant (rs1799971) upregulates OPRM1 expression across reward networks. Other signals of interest included IL6R and SHISA9 (EUR); GHR (AFR); and ASTN2 (AMR). Conclusions: This study identifies DDX6 as a novel OUD risk locus, replicates associations with OPRM1 and FURIN, and highlights biologically plausible ancestry-specific signals in AFR and AMR populations. We also replicated top variants in an independent population. Finally, integrating GWAS with deep-learning annotations provides specific, localized biological hypotheses to guide future experimental validation and targeted therapeutics.

Lithium has been the gold standard for the treatment and prevention of relapse in bipolar disorder for over 60 years. Guidance from the National Institute for Health and Clinical Excellence states explicitly to 'offer lithium as a first-line, long-term pharmacological treatment for bipolar disorder'. Yet, in the last two decades its use has been in decline with clinicians favouring anticonvulsants or antipsychotics when treating this condition. In this study, we have used three openly available datasets containing prescribing data from primary and secondary care to explore trends in the use of lithium in England, showing both regional and temporal variance between 2015-2024. We have shown that lithium use declined in primary care by 20.9% in the last ten years (2015-2024) and 10.9% overall in the last five years (2019 to 2025). We have also shown how there is some regional variation in the source of lithium for patients, although the vast majority is prescribed in primary care. Further research into clinical behaviour is needed to understand what is driving the decrease in lithium usage, and what barriers and enablers may influence its use across the country.

Estimating dengue force of infection (FOI) is essential for understanding transmission dynamics and targeting intervention programmes, yet surveillance data in endemic settings required for estimations are often incomplete, with varying formats. We developed a Bayesian hierarchical catalytic model that jointly fits age-stratified case data, aggregate case data, and seroprevalence surveys within a single framework, incorporating external covariates to improve parameter identifiability. Synthetic validation showed that covariates alone recovered accurate FOI point estimates even when most districts contributed only aggregate data, but did so with poorly calibrated uncertainty; anchoring the model with a single seroprevalence survey was necessary to bring credible interval coverage close to nominal. Applied to 128 districts across Java and Bali, Indonesia (2016-2024), the model revealed substantial spatial heterogeneity in FOI and reporting rates. Many districts in Java exceeded the WHO-suggested seroprevalence threshold for vaccine introduction, yet were classified as low-priority when using reported incidence as prioritisation criterion, particularly in areas with weak surveillance. Model-based seroprevalence estimation, integrating multiple data sources, offers a more consistent basis for identifying high-priority districts for vaccine introduction, and is less susceptible to surveillance bias than reported incidence.

Background: Children born prematurely (before 37 weeks) or admitted to the neonatal unit (NNU) are at increased risk of adverse long-term physical health outcomes. It is also recognised that there is an association with later academic performance and special educational needs, however it is not clear whether these broad risk factors could be used as stand-alone heuristics to identify children who may benefit from additional support in educational settings. We aimed to examine the associations between neonatal unit (NNU) admission and educational attainment in mid-childhood. Methods and Findings: Pregnancy data from a prospective birth cohort (Pregnancy Outcome Prediction Study, Cambridge, United Kingdom, 2008-2012) were linked to national educational outcomes (Department for Education, United Kingdom). Multivariable regression models adjusted for maternal, child, and socioeconomic factors were used to evaluate associations between (i) all NNU admissions, (ii) at term NNU admissions >48 hours, (iii) preterm birth without ongoing physical health needs, and educational outcomes at ages 5-11 years. Children who required any NNU care were more likely not to meet expected educational standards across multiple ages and domains in early and mid-childhood: age 5 early year foundation (aOR 1.64, 95% CI 1.19-2.27, p=0.003), phonics at age 6 (aOR 2.43, 95% CI 1.72-3.57, p<0.001), and at age 7 (here assessments were divided into multiple domains): reading (aOR 1.67, 95% CI 1.18-2.38, p=0.004), writing (aOR 1.72, 95% CI 1.25-2.38, p<0.001), mathematics (aOR 1.56, 95% CI 1.09-2.22, p=0.020), and science (aOR 1.85, 95% CI 1.22-2.78, p=0.003). Similar patterns were observed among both at term-born infants who stayed >48hrs in NNU (phonics assessment at age 6 aOR 2.26, 95% CI 1.51-3.36, p<0.001) and in children born preterm without long-term physical health sequelae (phonics assessment at age 6 aOR 3.07, 95% CI 1.96-4.81, p<0.001). These associations were robust to adjustment for demographic, perinatal, and socio-economic factors. By age 11, differences in academic attainment were attenuated and no longer clearly distinguishable across all exposure groups. However, there was an increased likelihood of special educational needs (SEN) at age 11 associated with any NNU admission (aOR 1.78, 95% CI 1.15-2.73, p=0.009), at term NNU admission for >48hrs (aOR 1.88, 95% CI 1.19-3.00, p=0.007), and children born preterm without long-term physical health sequelae (aOR 1.50, 95% CI 1.00-2.25, p=0.049). Predictive performance of any NNU admission for SEN at age 11 was moderate (AUC 0.70, 95% CI: 1.14-2.65, p=0.010), with balanced sensitivity and specificity and high negative predictive value. Conclusions: NNU admission, for both term and preterm infants, is associated with poorer educational outcomes and an increased likelihood of special educational needs in mid-childhood.